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First-principles study of the Kondo physics of a single Pu impurity in a Th host 
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Based on its condensed-matter properties, crystal structure, and metallurgy, which includes a phase diagram 
with six allotropic phases, plutonium is one of the most complicated pure elements in its solid state. Its anoma¬ 
lous properties, which are indicative of a very strongly correlated state, are related to its special position in the 
periodic table, which is at the boundary between the light actinides that have itinerant 5/ electrons and the heavy 
actinides that have localized 5/ electrons. As a foundational study to probe the role of local electronic corre¬ 
lations in Pu, we use the local-density approximation together with a continuous-time quantum Monte Carlo 
simulation to investigate the electronic structure of a single Pu atom that is either substitutionally embedded in 
the bulk and or adsorbed on the surface of a Th host. This is a simpler case than the solid phases of Pu metal. 
For the Pu impurity atom we have found a Kondo resonance peak, which is an important signature of electronic 
correlations, in the local density of states around the Fermi energy. Furthermore, we show that the peak width 
of this resonance is narrower for Pu atoms at the surface of Th than for those in the bulk due to a weakened Pu 
5/-ligand hybridization at the surface. 

PACS numbers: 74.25.Jb, 74.20.Pq, 71.27,+a, 71.28.+d 


I. INTRODUCTION 

Pure plutonium is one of the most exotic elemental met¬ 
als with respect to its condensed matter properties, crystal 
structure, and metallurgy [1-6]. It has a phase diagram with 
six allotropic phases, with the low-temperature monoclinic a- 
phase stable up to 395 K, and the technologically important 
face-center-cubic (fee) (5-phase stable between 592 Kelvin and 
724 Kelvin. While Pu metal exhibits a significant atomic- 
volume expansion, with the (5-phase 25% larger in volume 
than the o-phase, no magnetism has been experimentally ob¬ 
served for the (5-phase volume [7], These unusual physical 
and mechanical properties [3] are believed to be caused by the 
itinerant-to-localized crossover of 5/ electronic orbitals [4], 
Band-structure calculations within the local density approx¬ 
imation (LDA) fail to predict the equilibrium volume of the 
nonmagnetic (5-phase of Pu [8, 9]; while the LDA +U calcu¬ 
lations gives the experimental (5-phase volume with an appro¬ 
priate value of U but also indicates an instability toward an 
antiferromagnetic ground state [10-12], Later LDA+t/ calcu¬ 
lations [13-16] that included spin-orbit coupling (SOC) sug¬ 
gested that the magnetism could be naturally quenched and 
that the correct 5-phase volume could be obtained if the Pu 
atom is in a closed-shell Pu 5 f 6 configuration,. However, the 
LDA+L' method is unable to obtain the major 5/-character 
quasiparticle peak near the Fermi energy as observed by pho¬ 
toemission spectroscopy (PES) on 5-Pu [17, 18]. Recently, 
the LDA+DMFT method [21-23], which is a combination of 
LDA with the dynamical mean-field theory (DMFT) [19, 20] 
needed to capture important quantum dynamical fluctuation 
effects, has become a powerful tool to address the strongly 
correlated electronic materials. When applied to the 5-phase 
of Pu, the LDA+DMFT method can explain not only the large 
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volume expansion [22, 24] but can also address the absence of 
magnetism [25, 26] and the quasiparticle peak near the Fermi 
energy [26-31]. The same approach has also recently demon¬ 
strated site-selective electronic correlation effects in the a- 
phase of Pu [32]. The key to this success lies in the fact that 
the DMFT generates both the incoherent as well as the coher¬ 
ent parts of the spectral density in the intermediate coupling 
regime, with the former responsible for atomic-like behav¬ 
ior while the latter for the quenching of magnetism through 
Kondo screening [33], 

Experimentally, PES has played a significant role in uncov¬ 
ering the correlated-electron phenomena in /-electron materi¬ 
als, including elemental Pu solids. This probe measures the 
occupied part (that is, that below the Fermi energy) of the 
spectral density. However, to understand the precise nature 
of the quasiparticle bands in 5-Pu and other Pu-based ma¬ 
terials, it would be important to map out the quasiparticle 
states in the unoccupied regions as well. Scanning tunnel¬ 
ing microscopy (STM) has proved to be another important ex¬ 
perimental technique for understanding emergent phenomena 
in strongly correlated electronic systems such as, for exam¬ 
ple, high-temperature superconductors [34, 35]. This probe 
measures the local tunneling conductance at the atomic scale. 
One of the earlier successes of this powerful technique was 
the imaging of the Kondo resonance state for magnetic adap¬ 
tors (like Co) on metal surfaces [36-38]. Recently, this tech¬ 
nique has also been successfully extended to 4/-based heavy 
fermion systems [39, 40] and URu 2 Si 2 [41, 42], opening up 
new avenues of approach towards an understanding of 5/- 
electron properties in Pu. No doubt, a clear-cut STM obser¬ 
vation of Kondo resonance around a Pu impurity in a non¬ 
magnetic host would provide direct evidence for a correlation- 
driven entanglement of the Pu 5/-electrons with noninteract¬ 
ing conduction electrons in a metallic environment. 

While much theoretical work has been done on the elec¬ 
tronic correlation effects of the solid phases of Pu, the Kondo 
physics of a single Pu atom that is either embedded substi- 
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FIG. 1. (Color online) Geometries of a Pu impurity substituted in a 
Th bulk host (left panel) and absorbed on a Th (001) surface (right 
panel). The Pu atom is marked as a red ball while those of Th are 
blue balls. The XCrySDen program [48] was used to generate this 
figure. 


tutionally in the bulk or adsorbed on the surface of a tho¬ 
rium (Th) host is actually very important as a foundational 
study for Pu correlation physics, since such a study provides a 
benchmark for the electronic properties of an isolated Pu im¬ 
purity in an otherwise non-interacting bath. In addition, an 
impurity-like physics plays a critical role within the DMFT 
methodology, which makes a study of an isolated Pu impurity 
atom even more valuable. With this in mind, in this article 
we present calculations of the local electronic structure of Pu 
impurities in both the bulk and absorbed on the surface. This 
quantum many-body problem is treated by a combination of 
the local-density approximation with a continuous-time quan¬ 
tum Monte Carlo method. By considering the competition 
between the hybridization of the 5/ electrons with the de¬ 
localized conduction electrons and the local Coulomb repul¬ 
sion within the Pu-5/ electron manifold, Kondo resonances 
are found in the local Pu-5/ spectral density. By comparing 
the width of the resonance peak between the geometries, we 
find the Kondo temperature is higher for Pu atoms substitu- 
tionally embedded in the Th bulk than for those adsorbed on 
the surface of the Th host. Also, due to the strong spin-orbit 
coupling in the system, the resonance occurs in the j = 5/2 
manifold. 


From the perspective of electronic and magnetic properties, 
the Pu impurity problem is especially important because of its 
flexibility to address the electron localization/delocalization 
of Pu 5/ electrons. In addition, such an understanding will 
not only serve as a test of sophisticated theoretical capabili¬ 
ties but will also stimulate local-probe experiments like STM 
to contribute insights into emergent phenomena in Pu-based 
materials. 


II. THEORETICAL MODEL AND COMPUTATIONAL 
METHODS 

A realistic description of the electronic properties of a Pu 
atom embedded in the bulk or adsorbed on the (001) surface 
of a Th host, with the geometry shown in Fig. 1, can be for¬ 
mulated by a multi-orbital Anderson impurity model, with 

H = Hq + 'Hhyb + 'Hloc (1) 

and 

H 0 = X] e kd{d k , (2a) 

k 

^hyb — J2( v k a dlf a + h.c.) , (2b) 

k,a 

H\oc = '}2e a fif a + ^ HccP'ysflfpfsf'i ■ (2c) 

a; af3-y8 

Here, d\ (d k ) are the creation (annihilation) operators of the 
conduction-electron continuum for the state k of the Th bath 
while fl (fa) are those operators of the Pu impurity / elec¬ 
trons with quantum label a. The quantity e k is the energy dis¬ 
persion for the Th bath, where the label k includes crystal mo¬ 
mentum, band index, and spin, e a is the on-site single-particle 
energy levels, where a denotes both the orbital and spin quan¬ 
tum numbers, U a p-y$ are the Coulomb matrix elements, and 
14 a is the strength of hybridization between the local degrees 
of freedom and the continuum. 

Except for the local Coulomb interaction, our density- 
functional-theory (DFT) calculations for the Anderson im¬ 
purity model Hamiltonian (see just above) were done in the 
generalized gradient approximation (GGA) [43] using the 
full-potential linearized augmented plane wave (FP-LAPW) 
method of the WIEN2k code [44]. The spin-orbit coupling 
was included in a second variational step, with relativistic 
P 1/2 local orbitals added to the basis set for the 6 p states [45] 
of plutonium and thorium. The muffin-tin radius was set to 
3.3ao for Pu in the Th bulk and 2.5ao for Pu on the Th surface, 
while 2.5ao was used for all Th atoms in both cases. Here ao 
is the Bohr radius. The 5x5x5 number of k-points was 
used for the calculations of Pu in the Th bulk while 7x7x3 
for the calculations of Pu on the Th surface. We constructed 
the supercell by starting with a 2 x 2 x 2 number of fcc-Th 
conventional cells with the experimentally determined lattice 
constant 5.12 A [46]. For the Pu on the Th bulk, Pu was 
substituted for one Th (see the left panel of Fig. 1) while for 
the Pu on the Th surface, Pu was optimized along the z direc¬ 
tion toward the center of a Th surface plaquette (see the right 
panel of Fig. 1). The optimal distance of Pu away from the 
Th surface is 1.85 A in a vacuum of a 10 A perpendicu¬ 
lar distance between two neighboring Th slabs. We note that, 
because the calculation of systems with a large number of ac¬ 
tinide atoms is computationally demanding in the FP-LAPW 
method, a full-structure relaxation procedure is not explored. 
Furthermore, since the atomic number of Th is close to that 
of Pu, we expect that a full-structure relaxation will not qual¬ 
itatively change the Kondo physics discussed in the present 
work. 






3 


To include in the DFT an explicit of Coulomb interaction 
between Pu 5/-electrons requires a clear definition of the 
atomic-like local orbitals. In the present work, we use the 
weight-conserved projection procedure [47] to extract the lo¬ 
cal Green’s function for the correlated Pu 5/-orbitals from the 
full Green’s function defined in the DFT basis. Since only the 
5/ electrons of the Pu impurity site have an onsite Coulomb 
interaction, we can integrate out the non-interaction bath de¬ 
grees of freedom from the Th media and arrive at an effective 
action 

S e &=-^2[ [ drdr' /a(r)^l 0 " / l t^/ (r,r , )/ 1 ,(r , ) 

fny J o do 

ft* 

+ / drUiocif^T), f v (r)] . (3) 

J o 

Here, the fermion operators in the Hamiltonian become Grass- 
man variables in the action, and /3 = 1/T is the inverse of the 
temperature T. The non-interacting impurity Green’s function 
is 

= (ioJn A^j,(icU n ) (4) 

with a hybridization function 

a a. . \ _ N~^ VktiYkv 

— y , ( 5 ) 

“ *w„ - Ck 

where we use a Matsubara frequency formalism, with uj n = 
(2 n + 1 )T (n = 0,1,...) for fermions. A strong-coupling 
version of the continuous-time quantum Monte Carlo (CT- 
QMC) method [49-51], which provides numerically exact so¬ 
lutions, was used to solve this multiple-orbital Anderson im¬ 
purity problem. Although the Coulomb interaction parame¬ 
ters can also be determined by DFT calculations, we used a 
value of U = 4 eV, which is consistent with previous work 
on elemental Pu [24, 26, 29, 52]. The remaining Slater inte¬ 
grals, F 2 = 6.1 eV, F 4 = 4.1 eV, and F 6 = 3.0 eV, were 
calculated using Cowan’s atomic-structure code [53], but re¬ 
duced by 30% to account for screening. Since the DFT al¬ 
ready includes the Hartree-term of the Coulomb interaction, a 
double-counting correction Ed c = U(n ( j — 1/2) — J(n° —1)/2 
with a nominal value of = 5 for Pu-5/ electrons [52] has 
been used in the single-particle /-level. For very-well de¬ 
fined local Pu-5/ orbitals, this double-counting scheme has 
the virtue of numerical stability [54]. The Pu-5/ spectral den¬ 
sity from the CT-QMC simulations is extracted using the max¬ 
imal entropy method [55]. A similar procedure has been used 
to study the problem of a transition-metal magnetic atom Co 
in Cu hosts [56], where the formation of a Kondo-resonance 
Fermi liquid is uncovered. In that kind of systems, since the 
spin-orbit coupling is negligible, the cubic harmonics are a 
good basis to express the self-energy and local Green’s func¬ 
tion of correlated Co-3 d orbitals. In the present work, due to 
the strong spin-orbit coupling of actinide atoms, we choose 
the relativistic harmonics jrrij basis to express the self-energy 
and local Green’s function of correlated Pu-5/ orbitals. This 
choice minimizes the off-diagonal elements of the self-energy 
and local Green’s function, which reduces the minus-sign 
problem in QMC impurity simulations. 



FIG. 2. (Color online) DFT-based local 5/ projected density of states 
(DOS) of the Pu impurity as well as the 7 p, 6 d, and 5/ projected 
DOS of the Th sites nearest-neighboring to the Pu impurity for the 
cases (a) of Pu in bulk Th and (b) of Pu on the Th surface. The 
vertical dashed line in the figure denotes the Fermi energy. 


HI. RESULTS AND DISCUSSION 

In Fig. 2, we show the projected density of states (DOS) on 
the Pu impurity and that on Th nearest-neighbor sites around 
the impurity, obtained from DFT-GGA calculations. For the 
Pu atom in bulk Th, the strong spin-orbit coupling of Pu 
causes the 5/ states to be split into two manifolds, correspond¬ 
ing to j = 5/2 and j = 7/2 of the total angular-momentum 
quantum number. These states are located around the Fermi 
energy in the range from about -1 eV to 0.5 eV for the j = 5/2 
states and from about 0.5 eV to 1.5 eV for the y = 7/2 states. 
The Th-5/ partial DOS also exhibits SOC-split peak structure, 
but about 2 eV above the Fermi energy, suggesting that the Th- 
5/ states are not occupied. The Th-6d and Th-7p DOS show 
wide-band behavior, while the Th-6<7 contributes dominantly 
to the bath DOS around the Fermi energy, which hybridizes 
most strongly with the Pu-5/ impurity electronic states. For 
the Pu impurity on the Th surface, qualitatively similar fea¬ 
tures are exhibited in the DOS. Noticeably, the Pu-5/ DOS 
shows a narrower quasiparticle peak width with the dominant 
spectral weight distributed in the range of about -0.6 eV to 0.1 
eV for they = 5/2 states and 0.5 eV to 1.2 eV for they = 7/2 
states. This difference is closely related to the unique struc¬ 
ture of each geometry. For substitutional Pu in bulk Th, the Pu 
atom has 12 nearest-neighbor Th sites, with a Pu-Th distance 
of 3.60 A; while for Pu atoms on the Th surface, there are 
only 4 bonds to the nearest Th atoms, with a bond distance of 
about 3.15 A. 

To further clarify the close relation between the electronic 
structure and the local atomic environment for each geometry, 
we evaluate the frequency-dependent hybridization function, 
r(w) = — Im[A(w)]/7r, where A(w) is the hybridization self¬ 
energy given by Eq. (5). We show in Fig. 3 the hybridization 
function obtained from the GGA calculations. For the Pu im- 
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FIG. 3. (Color online) Hybridization function between the Pu-5/ 
states and the Th bath calculated within the GGA approximation for 
(a) Pu in bulk Th and (b) for Pu on the Th surface. The vertical 
dashed line in the figure denotes the Fermi energy. 


purity in bulk Th, the hybridization function strength is about 
0.8 eV when averaged over the frequency range between -0.3 
eV and 0.7 eV. For the Pu atom on the Th surface, the strength 
is about 0.5 eV when averaged over the frequency range be¬ 
tween -1 eV and 1 eV. In addition, the overall hybridization 
function strength for both the j = 5/2 and the j = 7/2 man¬ 
ifolds is similar, with that for j = 7/2 being only slightly 
larger for a given geometry. We ascribe this slight difference 
to the fact that the Pu-5/ level for the j = 7/2 states, e J=7 / 2 , 
is closer to the center of gravity of the Th-(k/ band. 

By switching on the Coulomb interaction on the Pu impu- 
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FIG. 4. (Color online) Total angular-momentum-resolved density of 
states of Pu impurity in Th (a) and of Pu on the Th surface (b) ob¬ 
tained from the CT-QMC simulations at temperature T = 232 Kelvin. 


rity, the CT-QMC simulations indicate a significant electronic 
correlation on the Pu impurity. For the Pu atom in bulk Th, 
as shown in Fig. 4(a), the total DOS (red line) obtained from 
the CT-QMC simulations at T = 232 Kelvin show not only 
a quasiparticle peak very close to the Fermi energy but also 
incoherent peaks at higher energies. From the plot of the j- 
resolved DOS, one can see clearly that the j = 5/2 channel 
(blue line) is the dominant contribution to the quasiparticle 
peak with only minor intensity from the j = 7/2 channel at 
the same energy position (green line), while the incoherent 
peak located above the Fermi energy appears mainly in the 
j = 7/2 channel (green line). The quasiparticle peak obtained 
from the CT-QMC simulations is also more than a factor of 2 
narrower than the peak on the j =5/2 partial DOS obtained 
from GGA calculations, suggesting a many-body interaction- 
driven Kondo resonance state. For the Pu atom on the Th 
surface, which has a reduced hybridization, the quasiparticle 
peak becomes much narrower, as shown by the red and blue 
lines in Fig. 4(b), especially when compared with the case 
of a substitutional Pu atom in bulk Th. In addition, the in¬ 
coherent peak below the Fermi energy also becomes visible 
in the j = 5/2 channel. Therefore, our results indicate that 
the Kondo resonance state develops in the j = 5/2 channel. 
Physically, it results from almost all of the occupied Pu-5/ 
electrons residing in the j = 5/2 subshell. When compar¬ 
ing the results from GGA calculations (cf. Fig. 2) and from 
the low-temperature CT-QMC simulations (cf. Fig. 4), the 
degree of interaction-driven quasiparticle renormalization be¬ 
tween the two geometries is similar. To quantify the quasipar¬ 
ticle renormalization, we estimate the quasiparticle spectral 
weight defined as 


1 i9ImS(*u; rl ) 


dbJn 

w n =0- 


where the small contribution from the hybridization at the 
Fermi energy has been neglected. Because the quasiparticle 
peak occurs in the j = 5/2 subshell, we consider only the 
correlation-driven self-energy in this channel, which leads to 
Zj— 5/2 ~ 0.33 for the Pu in bulk Th while Z ss 0.35 for the 
Pu on the Th surface. When the numerical uncertainty is taken 
into account, the above observation is well supported. 

Within the CT-QMC simulations, the Kondo temperature 
that is defined from a renormalized fermi-liquid theory for the 
Anderson impurity model [57] that is given by 

T k = - 7 rZImA( 0)/4 (7) 

can also be determined. It was found that Tk ~ 968 Kelvin 
for the Pu in bulk Th and T K ss 457 Kelvin for the Pu on 
the Th surface. This energy scale should be scaled with the 
width of the Kondo resonance peak developed near the Fermi 
energy, which can be measured by STM experiments. To 
demonstrate the temperature effect, we also performed CT- 
QMC simulations at a higher temperature (T = 725 Kelvin). 
As shown in Fig. 5, the Pu-5/ spectral density for the case 
of the Pu in bulk Th is not much different than that at low 
temperature while the quasiparticle peak for the case of the 
Pu on the Th surface is strongly suppressed in intensity. In¬ 
terestingly, at this temperature, which is already higher than 
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FIG. 5. (Color online) The same as in Fig. 4 except at the temperature 
T = 725 Kelvin. 

the Kondo temperature for the Pu on the Th surface, the low- 
energy DOS peak is still quite visible. One possibility is 
that the Kondo temperature is systematically underestimated. 
The other possibility is that the low-energy DOS peak may 
contain the contribution from both the Kondo-driven many- 
body physics and the multiplet effect usually occurring in the 
atomic physics. Our CT-QMC simulations yield a Pu-5/ oc¬ 
cupancy of (n/) ss 5.29 for the Pu in bulk Th and 5.35 for 
the Pu on the Th surface. These values suggest that the Pu-5/ 
electronic states are in the mixed-valence regime. The Pu- 
5/ mixed-valence behavior together with the multiplet effects 
have also been proposed in the elemental Pu solid [26, 27]. 
Therefore, this resemblance leads one to believe that the sec¬ 
ond possibility is more likely. Finally, it is worthwhile to men¬ 
tion that we also performed the CT-QMC simulations for the 
case of a Pu impurity substituted for a Th atom in the sur¬ 


face of the Th slab. In this geometry, the Pu impurity has 8 
nearest-neighboring sites with the Pu-Th distance of 3.60 A. 
We obtained results (not shown here) similar to the case of the 
Pu on the Th surface. 


IV. SUMMARY 

We have performed a first-principles study of the Kondo 
physics of a Pu impurity in a Th host within the framework 
of DFT theory in combination with the CT-QMC simulations. 
By considering two representative geometries of either a sub¬ 
stitutional Pu impurity embedded in bulk Th bulk or adsorbed 
on the Th surface, we have been able to identified the Kondo 
resonance states developing near the Fermi energy. For a fixed 
set of on-site Coulomb-interaction parameters, the width of 
the quasiparticle DOS peak has been shown to be sensitive 
to the strength of the low-energy hybridization of Pu-5/ elec¬ 
trons with the conduction band in the Th medium, as borne out 
by the two geometries considered here. By varying the tem¬ 
perature, we have also discussed the nature of the quasipar¬ 
ticle peak in the context of the estimated Kondo temperature 
scale. With the recent success of STM for exploring emergent 
phenomena in /-electron materials, we hope that the results 
presented in this work will stimulate STM exploration for Pu- 
based /-electron materials. 
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